Low-dimensional Bose liquids: beyond the Gross-Pitaevskii approximation 
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The Gross-Pitaevskii approximation is a long- wavelength theory widely used to describe a variety 
of properties of dilute Bose condensates, in particular trapped alkali gases. We point out that 
for short-ranged repulsive interactions this theory fails in dimensions d < 2, and we propose the 
appropriate low-dimensional modifications. For d — 1 we analyze density profiles in confining 
potentials, superfluid properties, solitons, and self-similar solutions. 
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Experimental observation of Bose-Einstein condensa- 
tion (BEC) in trapped alkali vapors [l] has ushered in 
a new era of superlow temperature physics. The Gross- 
Pitaevskii (GP) mean-field theory || has proven to be an 
indispensable tool both in analyzing and predicting the 
outcome of experiments. 

With rapid progress in the experimental exploration of 
BEC systems it is reasonable to anticipate that effectively 
one and two-dimensional systems are realistic prospects 
in the near future ||. For example, high aspect ratio 
cigar-shaped traps approximating quasi-one-dimensional 
BEC systems are already available experimentally |Q. 
From the theoretical viewpoint low-dimensional systems 
are rarely well represented by mean-field theories which 
leads one to question the validity of the GP theory in 
one and two dimensions. In this paper we shall show 
that, indeed, the physics of dilute Bose systems requires 
a fundamental modification of the GP theory in low di- 
mensions d < 2. 

The GP theory is a quasi-classical (or mean-field) ap- 
proximation which replaces the bosonic field operator ip 
by a classical order parameter field $(r, t). For short- 
ranged interactions the interparticle potential U(r) is re- 
placed by g5 d (r), where g is the pseudopotential. Then 
for a system of bosons in an external potential V, the 
energy functional and the equation of motion for $ take 
the following form: 
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The GP equations (|l|) and (||) are widely used to compute 
a variety of properties of Bose systems . 

The GP approximation is a long-wavelength theory re- 
lying on the concept of the pseudopotential to account for 
interparticle interactions. However, for repulsive bosons 
the canonical pseudopotential vanishes in two dimensions 
||, implying that an essential modification of the GP 
theory is necessary for d < 2. To see how to modify 



the theory, it is useful to rewrite the last integrand of 
(|l|) in terms of the particle density n = |$ 2 , so that 
(,g/2)|$| 4 = {g/2)n 2 . This can be recognized §,§ as the 
lowest order term of a dilute expansion of the ground- 
state energy density for d > 2. Thus the correct low- 
dimensional local density theory will instead have the 
ground-state energy density of the d < 2 dilute Bose sys- 
tem HH, which is not proportional to n 2 . 

Let us write the interparticle interaction in the form 
U(r) = u S d (r) where uo is the amplitude of the inter- 
particle repulsion, and the notation S d (r) denotes any 
well-localized function that transforms into the mathe- 
matical Dirac S function when the range of interaction 
a — > 0. The renormalization group analysis of this prob- 
lem B reveals that two dimensionless combinations na d 
and h rS 2 ~ d ^ d /muQ play an important role in defining 
the conditions of the dilute limit for d < 2. 

In the dilute limit na -C 1 and h 2 n/muQ -C 1, any 
one-dimensional Bose system with short-ranged repul- 
sive interactions becomes equivalent to a gas of free- 
fermions @ (or equivalently point hard-core bosons ||), 
with energy density 7r 2 ?i 2 n 3 /6m. This can be general- 
ized [[| for d < 2 to the statement that for na d <C 1 
and h n( 2 ~ d ^ d /muo <C 1, the lowest order term of the 
ground-state energy density expansion is universal and 
given by (h 2 Cd/2m)n ( - 2+d ' ) / d , where Cd is a d-dependent 
constant. This implies that the quartic nonhnearity |<I>| 4 
in Eq.(||) should be replaced by |<J>| 2 ( 2 + rf )/ d . In particu- 
lar, for the practically important case of one dimension, 
the system of equations (|lj) and ^ is modified to 
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Similar considerations applied to the marginal two- 
dimensional case lead to the conclusion that in the di- 
lute limit |^,|6| a theory replacing the GP approximation 
starts from the energy functional 
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Ignoring the logarithmic factor will perhaps suffice for 
many practical purposes; then (^|) is precisely of the GP 
form. Finding an experimental manifestation of the log- 
arithmic correction is likely to be very challenging. 

A detailed analytical and numerical study of the one 
and two-dimensional cases will be given in a longer pub- 
lication ||; hereafter we restrict ourselves to only the 
salient features of one dimension where the deviations 
from the GP theory are largest. 

Density profiles in external potentials: The station- 
ary solution to Eq.(||) defined via &{x,t) — (ft(x)e~ tpit / h 
can be found by solving 
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subject to the condition of fixed total particle number 
N = J dx(f> 2 which determines the chemical potential \i. 
For an external potential that varies slowly on the scale 
of the interparticle spacing the derivative term in (||) can 
be ignored: this gives the density profile in the Thomas- 
Fermi (TF) approximation: 
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(x) = ct> 2 TF = [2m( f i - V(x))} 1/2 /nh 
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with the density being zero in the classically forbidden 
region /i < V(x). For the practically important case of a 
harmonic trap, V — muj 2 x 2 /2, and the density profile is 
elliptical: 
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The chemical potential is given by [Itf = hojN, 
the density in the center of the trap is ?t-tf(0) = 
{2mhu)N) 1 / 2 /nh, and the size of the trapped condensate 
is 2(2mhojN) 1 / 2 /irh. 

The accuracy of these predictions can be tested against 
the exact solution of a dilute system of bosons with re- 
pulsive interactions: an ideal candidate being a system 
of point impenetrable bosons. The boson-fermion equiv- 
alence Q] implies that in the many-body system, the 
single-particle energy levels E n — hui(n + 1/2) of the 
harmonic oscillator are occupied in a fermionic fashion, 
i.e. with no more than one particle per state. The chemi- 
cal potential is then given by fj, = hcu(N + 1/2), which for 
large N approaches our TF result (j,tf = hcuN. Similarly 
the density profile can be computed as a sum of squares 
of the single-particle wave functions: 
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where Hk are Hermite polynomials and I = (fi/muj) 1 / 2 . 
The density distribution (|9|) is plotted in Fig.l where it 
is compared with a) the numerical solution of (Q) with 
V = muj 2 x 2 /2 1 and b) the TF result (§), for different 
numbers of particles. The main flaw of the theory based 
on (|^) is that it does not reproduce density oscillations 
due to algebraic ordering of the particles. This is not 
surprising as (akin to the GP approximation) the dis- 
creteness of the particles, which is responsible for the 
density oscillations, is ignored. Otherwise, the agree- 
ment between the approximate and the exact profiles 
is very good; in the limit of large particle number the 
differences become imperceptible. These results can be 
directly tested experimentally; as a comparison we note 
that the one-dimensional GP theory in the TF approxi- 
mation predicts Q utf ~ — V(x), which is quite dis- 
tinct from (Q), and agrees very poorly with the exact 
result. 




a) b) 
FIG. 1. The density profile (^J) plotted for particle num- 
bers N = 4, 9, 16, 25 in units mui/h — 1. The non-oscillating 
curves correspond to a) the numerical solution of (^) and b) 
the TF result (|). 

Solitons: Gray solitons |l(J have been recently created 
and their dynamics was observed in cigar-shaped con- 
densates of 87 Kb vapors which makes it important 
to understand solitonic properties of the system (^|) and 
(|). Let us look for solutions to (|) (with V = 0) of the 
form <I>(.t,<) = (j>{x, t)e~ %lit l h . The function then obeys 
the equation 



ihd t 



2m 



(10) 



where the chemical potential /i = ■K 2 % 2 4'f ) /2m is selected 
so that the particle density uq = </>q is constant at infinity. 
In dimensionless variables / = 4>/4>o, y = imox, t = 
Tr 2 nlht/m, Eq.(|i~0|) simplifies to 



2id T f = -d 2 y f+(\f\ 4 -l)f 
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We will be looking for a localized traveling wave solu- 
tion |0 to (|ll|) of the form f(y,r) = f(y - (3t) where 
the dimensionless velocity (3 is measured in units of the 
sound velocity c = Trhno/m. This problem can be solved 
exactly. The results are conveniently described in terms 
of the amplitude A and phase 9 of the dimensionless order 
parameter / = Ae : 
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The spatial behavior given by ( |l2] ) is shown in Fig. 2 
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FIG. 2. The density yl 2 and the phase 
soliton (|l|) with = 1/2. 



for the moving 



The amplitude in ( |l2[ ) describes a moving depression 
(particle deficit) with the minimal value at the soliton 
center given by A 2 (0) = (1 + 3/3 2 ) 1/2 - 1. The soliton 
exists only for f3 < 1 (i.e. the soliton velocity cannot 
exceed the speed of sound); for (3 = 1 ([l2]) gives the 
uniform result A 2 = 1. On the other hand, for (3 = 
(i.e. a vortex, or dark soliton [jTof ) the minimal value of 
the amplitude at the soliton center drops to zero. 

The phase expressed in ( |l2|) varies rapidly in the vicin- 
ity of the amplitude dip, staying approximately constant 
far away from it. The total phase shift across the soliton 
can be found as AO = co S - 1 [(3/3 2 - 1)/(1 + 3/3 2 ) 1 / 2 ]. It 
is a continuous function of the soliton velocity varying 
between ir (when j3 = 0) and zero (when (3=1). Anti- 
solitons may be defined as having opposite signs of dO/dy, 
and there are no constraints on AO for the open line or 
ring geometries (provided the number of solitons matches 
the number of antisolitons). However, if there is an im- 
balance of solitons and antisolitons in the ring geometry, 
then the uniqueness of the order parameter /(y, r) im- 
plies that AO is a fraction of 2ir for any excess soliton; 
this will in turn mean that the excess soliton velocity is 
quantized. 



The solution ( |12[) bears some similarity with the one- 
dimensional soliton of the GP theory Jl(| ; the main qual- 
itative difference (seen after recovering the original units) 
is that in the dilute limit the soliton size is of order 
1/(1 — /3 2 ) 1 / 2 rio independent of the amplitude of the in- 
terparticle repulsion Jh3[ . 

General methods |lCf| can be used to compute the soli- 
ton energy E and momentum P. For their dimensionless 
counterparts e = 2mE /i: 2 h 2 n^ and V = P/ttTluq we find 
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The dependencies e((3) and V(j3) parametrically define 
the soliton dispersion law e(V) which should be identi- 
fied |Q with the "hole" branch of the elementary ex- 
citations spectrum jl{|. To assess the accuracy of e(V) 
given in ( |l3| ) we compare it with the exact result of Licb 
p"5| for the system of ^-interacting bosons in the dilute 
limit h 2 n/mu < 1: e e xactCP) = 2\P\ - V 2 , for \V\ < 1. 
Since the velocity j3 in ([l3]) varies between zero and unity, 
the momentum (which we choose to be positive) com- 
puted from ( |l3| ) varies between unity and zero in corre- 
spondence with the exact result. It is straightforward to 
show that for P< 1, the elimination of (3 in (Jl^) leads to 
e = 2V which is again in agreement with the exact result. 
The behavior e(V) implied by ([l3|) in the vicinity of the 
end-point of the spectrum V = 1 is qualitatively similar, 
and quantitatively close to the exact dependence. To il- 
lustrate these statements we have plotted the dispersion 
law ( |l3| ) in Fig. 3 against the exact result. 




0.8 p 1.0 

FIG. 3. The spectrum parameterized by ( |l3| ) (lower curve) 
compared to the exact result of Lieb (upper curve). 

Super-flow: The dimensionless current density is given 
by j = A 2 d y 0, and below we look for solutions with fixed 
given current j (i.e. the steady state) and d T A = d T = 0. 
Substituting / = Ae w and j = A 2 d0/dy into (0), and 
imposing fixed chemical potential, we find: 
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In the spatially uniform state ^pr = and one finds the 
dimensionless amplitude A^ = [1 + (1 — 4j 2 ) 1 / 2 ]/2, which 
implies that superflow reduces the amplitude of the order 
parameter. The uniform solution, and thus superfluidity, 
cease to exist above the critical flow j c = 1/2 when the 
amplitude drops down to its minimal value A^ = 2~ 1 / 4 . 
These results imply that the critical velocity for super- 
fluidity in the original units is c/^/2. 

The equation ( |l4] ) also has an immobile well-localized 
solution in the form of a dip of the order parameter; far 
away from the dip the amplitude recovers to its uniform 



3 



value. The dip solution is closely related to the soliton 
previously discussed. Indeed, in the reference frame mov- 
ing with the flow, the dip solution is moving and thus is 
identical to a soliton. The functional form of the dip 
can be deduced from ( |l2| ) by replacing (3 by j/A^, A by 
A/Aao, and (y — (3t) by yA 2 ^. The dip solution disap- 
pears altogether for j > j c . 

Self-similar solutions: The results derived so far have 
their counterparts in the context of the one-dimensional 
GP approximation. However the theory based on Eqs.(||) 
and (^) allows self-similar solutions which do not exist in 
the one-dimensional GP theory |l6) . Below we only look 
at the cases consistent with the condition of conservation 
of total particle number. Consider the system of bosons 
placed in a harmonic trap. In dimcnsionless variables it 
is described by 
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where w — muj /■n 2 Tin 2 ) is the dimensionless oscillator fre- 
quency (no now has the meaning of a density introduced 
to make / dimensionless; it should be determined from 
the complete solution of the problem). In contrast to 
( pd| ) (and without loss of generality) we have shifted the 
origin of the chemical potential. The self-similar solution 
/ = Ae ie derived from (|l|) has the form jl6| 



A = p{rr 1 / 2 h(y/p{r)) 
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where the functions p(r) and h(v) obey the equations 
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where 7 and S are arbitrary constants. Eq.([l8]) for the 
scaling function h(v) has localized solutions only for i5 > 
and 7 > 0: for 7 = an explicit analytic solution to 
( |l8| ) can be written down, while for 7 > 0, ([l8]) has the 
same functional form as the equation we encountered in 
determining the density profile in the harmonic trap (cf. 
(|) for V = muj 2 x 2 /2}. 

The dynamics of the length scale p(r) can be under- 
stood by viewing (Jit]) as a fictitious classical mechanics 
problem in the potential U = w 2 p 2 /2 + j/2p 2 . This anal- 
ogy implies that an initially localized cloud of bosons in 
free space (w — 0) will expand asymptotically in a bal- 
listic fashion: p{r) ~ r. In the presence of the confining 
potential (w ^ 0) the scale p(r) oscillates between max- 
imum and minimum values: for 7 = the dynamics of p 
is the same as for a harmonic oscillator of frequency w. 

We have also performed direct numerical integration 
of the non- linear equation (||), and have confirmed the 
existence of both the similarity solutions, and moving 
trains of solitons with quantized velocity, with amplitude 



and phase as given by (|12|) . More details will be given in 
a future publication Q . 

In conclusion we have presented a new continuum de- 
scription of dilute Bose liquids appropriate for low di- 
mensional systems. For the case of one dimension, we 
have derived stationary properties, along with solitonic 
and similarity solutions. In particular, the latter have no 
analog in the GP theory. It is our hope that these results 
will be testable in BEC experiments in the near future. 
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